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Abstract 

We propose a flexible and model independent parametrization of the neutrino mixing 
matrix, which takes advantage of the fact that there are up to three small quantities in 
neutrino mixing phenomenology: (i) the deviation from maximal mixing of solar neutri- 
nos, (ii) the mixing matrix element C/gS and (iii) the deviation from maximal mixing of 
atmospheric neutrinos. It is possible to quantify those three observations with a parame- 
ter A ~ 0.2, which appears at least linearly in all elements of the mixing matrix. The limit 
A — > corresponds to exact bimaximal mixing. Present and future experiments can be 
used to pin down the power of A required to usefully describe the observed phenomenol- 
ogy. Observing that the ratio of the two measured mass squared differences is roughly 
A^ allows to further study the structure of the Majorana mass matrix. We comment on 
the implications of this parametrization for neutrinoless double beta decay and on the 
oscillation probabilities in long-baseline experiments. 



1 Introduction 



Neutrino physics has made impressive progress in recent years |lj. In particular, the structure 
of the neutrino mixing matrix has been identified to a reasonable precision. The final step 
for the determination of its structure has come from the KamLAND experiment |2I, which 
confirmed the large mixing angle solution for the solar neutrino problem, after a preference for 
this parameter space was already implied by the data of the various solar neutrino experiments 
in]. Very recently, the SNO salt phase data ^ rejected maximal solar mixing by more than 
5(j in El- The (almost) maximal mixing of atmospheric neutrinos has been found by the 
Super Kamiokande experiment [7] and confirmed by the K2K collaboration [S]. Finally, the 
presence of a small if not zero angle was implied by reactor experiments j^. In the present 
paper we wish to propose a parametrization of the neutrino mixing matrix in terms of a small 
parameter A, whose magnitude is interestingly around 0.2, i.e., close to the Wolfenstein pa- 
rameter used to parametrize the CKM matrix [1^. In any parametrization of the neutrino 
mixing matrix (for earlier attempts, see fllElEl), it is convenient to start from a reference 
matrix and describe deviations from it. Our reference matrix is the one corresponding to ex- 
act bimaximal neutrino mixing. The parameter A describes the up to three small deviations 
from this mixing scheme, namely the deviation from maximal mixing of solar neutrinos, the 
deviation from zero t/gs and the deviation from maximal mixing of atmospheric neutrinos. The 
magnitude of A is defined by the observed non-maximality of solar neutrino mixing |31 |H1 El 
and future precision experiments can be used to pin down the power of A to usefully describe 
the other two deviations. In addition, the ratio of the mass squared differences governing solar 
and atmospheric neutrino oscillations is given by A^, so that it is possible to analyze also the 
structure of the neutrino mass matrix (provided neutrinos are Majorana particles). We also 
analyze the oscillation probabilities for long baseline experiments and the effective mass as 
measured in neutrinoless double beta decay within our parametrization. 

The paper is organized as follows: In Section |21 we describe the neutrino mixing parameters as 
implied by current data and outline the idea of our parametrization. Then, in Section El we give 
the form of the mixing matrix for various special cases of the parametrization and analyze in 
Section 01 the form of the neutrino mass matrix. In Section |31 we apply our parametrization to 
the effective mass as measured in neutrinoless double beta decay and to long-baseline oscillation 
experiments. We conclude in Section (HI 

2 Quark versus Lepton Mixing 

The Wolfenstein parametrization JUl of the CKM matrix uses the fact that the quark mixing 
is very small, i.e., the mixing matrix is quasi the unit matrix with only small corrections on 
the off-diagonal entries. In terms of mixing angles, a hierarchy of the form 6*12 ~ 0.1 > 6'23 ~ 
0.01 > 6'i3 ~ 0.001 is observed. This has been used by Wolfenstein to introduce an expansion 
parameter A describing the mixing between u and s quarks. The observation that c — b {u — b) 
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mixing is roughly one (two) orders of magnitude suppressed then leads to 

/ 1-|A2 A AX^p-ir])\ 



CKM 



-A 



1-iA^ 



\AX^l-p + i7]) -AX^ 



AX^ 
1 



+ O(A^) 



Of course, A corresponds to the Cabibbo angle sin^^c" — 0.22, whereas the other parameters are 
about 14 A ^ 0.83, p ~ 0.23 and rj ~ 0.36. The latter parameter describes the CP violation 
in the quark sector; all such effects have to be proportional to 



JcP = lm{Vud V:, V:^} - -A^ A% ~ -3 ■ 10"^ 
Therefore, CP violation in the quark sector is a small effect. 



2.1 Neutrino Mixing 

The neutrino oscillation data can consistently be described within a 3-neutrino mixing scheme 
with massive neutrinos, in which the flavor states [a = e, p, r) are mixed with the mass 
states Ui {i = 1,2,3) via f/pMNS, the unitary Pontecorvo-Maki-Nagakawa-Sakata [Hj lepton 
mixing matrix. It can be parametrized as 



/ C12C13 S12C13 Si3 \ 

f^PMNS = "~'Sl2C23 — Ci2S23Si3e*'^ C12C23 — Si2S23Si3e*'' 5230136*'^ 

\ S12S23 - Ci2C23Si3e*'^ -C12S23 - Si2C23Si3e*'^ 0230136*'' / 



diag(l, e*°, e*^) , (3) 



where 5 is a Dirac CP violating phase, a and (3 are possible Majorana CP violating phases |17j . 
Cij = cos 6ij and Sij = sin 6ij. We shall not consider the two Majorana phases in this Section. 
The angles 612 and ^23 control the oscillations of solar and atmospheric neutrinos, respectively. 
The angle 613 is mainly limited by reactor z/g experiments: The Dirac phase S can be measured 
in long baseline neutrino oscillation experiments (see, e.g., |IHI)- 

To obtain information about the PMNS matrix one fits |51 El the results of neutrino ex- 
periments to the hypothesis of neutrino oscillations. The relevant formula for the oscillation 
probabilities is 

P{u^^up)=5^p-2ReJ2Ua^K^U;,Up, h-exp^^j , (4) 

j>i \ ^ 

where Am|j = — mf. The 1 (3)cr ranges of recent analyzes read in terms of the para- 
met riz at ion Eq. (jSjl: 

(0.27) 0.35 < tan2 Ou < 0.52 (0.72) |3] , 

(0.45) 0.75 < tan^ ^23 < 1-3 (2.3) , ^ (5) 
(0) < sin^ < 0.029 (0.074) . 
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The best-fit points are given by 

tan^^is = 0.43 0, tan^ ^23 = 1 [El • (6) 

For 6*23 there is an ambiguity corresponding to 6*23 ^ 7r/2 — 6*23, i.e., the angle hes on the 
"hght" or "dark side". Matter effects in future long-basehne experiments will resolve this. In 
terms of the often used parameter sin^ 26'23, which is blind to this ambiguity, one has at 1 (3) a: 
sin2 2e23 > 0.86 (0.84). 

Regarding the mass squared differences, the best-fit values are (AmQ)BF = 7.2 • 10~^ eV^ jH] 
and (Am^)BF = 2.6 ■ 10~^ eV^ ^Hl- A recent preliminary analysis of the SuperKamiokande 
collaboration, taking into account e.g. improved identification criteria and neutrino fluxes, yields 
a value of (Am^)BF = 2.0 ■ 10~^ eV^ There are two possible mass orderings, the normal 
and the inverse mass ordering: 

normal mass ordering: AttIq = Am2i ^ Am32 — Amg^^ = Am\ 

(7) 

inverse mass ordering: Ami, = Am2i ^ —Am'^i ~ — Am|2 = Am^ 

Extreme cases are the normal (inverse) hierarchy with = Am\ ^ m2 = Ami, ^ mi (m2 = 
Am^ ~ mi 3> ma) and the quasi-degenerate mass scheme for which m"^ ~ mg — mi ^ Am^. 
The latter is fulfilled for values of the neutrino masses larger than ~ 0.2 eV. 

Ignoring the phases, the "best-fit PMNS matrix" reads: 

/ 0.84 0.55 \ 

[/BF^g = -0.39 0.59 0.71 

\ 0.39 -0.59 0.71 / 
In the pre-SNO-salt-phase analysis of Ref. there was given the 3a range of the PMNS 



(9) 



/ 0.73 -0.88 0.47- 0.67 -0.23 \ 



If/i 



PMNS I 



0.17- 0.57 0.37- 0.73 0.56 - 0.84 



V 0.20 - 0.58 0.40 - 0.75 0.54- 0.82 / 
where the phase S was allowed to take arbitrary values. 

Next generation long-baseline experiments will be able to probe Am^ and sin 2^23 to the % 
accuracy The element Ue3 can be probed down to the level 10~^ in future long-baseline 
or reactor experiments j22]. Neutrino factories [22] can improve these bounds considerably. 
The solar neutrino mixing angle tan^ 6*301 will see its error reduced below 10 % by experiments 
investigating the low energy neutrino fluxes from the sun 

Currently no information about leptonic CP violation exists. In oscillation experiments one 
can detect CP violating effects which have to be proportional to 

Jcp = lm{Uei U^2 U*^ 11*2} = I sin 26*12 sin 26*23 sin 26*13 cos 6*13 sin 5 



Q 

(1 — corrections from 6*13 and non-maximal 6*12 and ^23) 



(10) 
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where the value ^13/4 is the hmit for small ^13, maximal 612 and ^23- 



2.2 The Strategy 

As a very useful limit, the bimaximal mixing pattern jSEl, corresponding to 612 = 623 = 7r/4 
and 6'i3 = 0, can be considered. The resulting mixing matrix, ignoring the CP violating phases, 
reads 



T rbimax 
•^PMNS 



/ 73 73 \ 

_ 1 1 1 

2 2 V2 

1 _ 1 J_ 

y 2 2 V2 y 



The non-zero entries therefore take values of and ±|. This form of f/pMNS shall be our ref- 
erence matrix, whose deviations are to be described by some small parameter A. In some sense, 
the matrix (lllj) is the analogue to the unit matrix in case of quark mixing. Corrections of order 
A and higher to the unit matrix lead to the CKM matrix as parametrized in Eq. (P). In the 
same way, corrections to f/pMNs order A and higher can lead to the observed neutrino mixing 
phenomenology with non-maximal solar neutrino mixing, possible non-maximal atmospheric 
neutrino mixing and non-zero f/es- One might state that the unit matrix in the quark sector 
and the bimaximal mixing matrix in Eq. (jllll are the zeroth order form of the relevant mixing 
matrix. We shall comment on a possible origin of corrections to bimaximal mixing (along the 
lines of Ref. fTf^) Section 13.71 Tries to parametrize the PMNS matrix in analogy to the 
CKM matrix suffer in general from the fact that from the 9 elements of ^7pmns only one is 
small, namely the element {Uesl ~ 0.27. The other eight entries take typically values around 
1/2 or I/V2. There are to our knowledge two other approaches to do something similar to the 
PMNS matrix as has been done so successfully with the CKM matrix. The analysis from ^I] 
uses Ue2 — \/2A, atmospheric neutrino mixing remains maximal and Ue3 is proportional to A^. 
The work ^] chooses the expansion parameter A = f/^3 ~ 1/ ^/2. Solar neutrino mixing is 
governed by Ue2 = AX^ and the small quantity f/es has to be introduced at at least the eights 
power of A. Effects due to CP violation in neutrino oscillations are — courtesy of Eq. (jiup — 
proportional to at least A^^. An Ansatz for t/pMNS corresponding to tan^^i2 = 0.5, maximal 
atmospheric mixing and zero Ues has been used in To describe deviations from it, it has 
been multiplied with a Wolfenstein-like matrix. 



In this letter we wish to propose a purely phenomenological and model independent para- 
metrization of the PMNS matrix by using a small "expansion" parameter A. For a useful 
analysis in terms of a small parameter one requires small quantities in the mixing matrix. The 
basic idea is given by the identification of up to three such small numbers in neutrino mixing 
phenomenology, namely: 

(i) the deviation from maximal mixing of solar neutrinos 

(ii) the small mixing element Ue^ 
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(iii) the possible deviation from maximal mixing of atmospheric neutrinos 

Those three aspects describe all possible deviations from the bimaximal mixing scheme in Eq. 
(fTTjl . Observation (i) is now a solid experimental evidence, after inclusion of the SNO salt phase 
data it now holds that tan^ ^0 < 1 at more than Sex [3101^]. Regarding observation (ii), 
only the mentioned limit of |?7e3p < 0.07 (at 3cr) exists. Best-fit points of three flavor analyzes 
of all neutrino data typically yield very small if not vanishing values for this quantity. Finally, 
atmospheric neutrino mixing insists to be described by solutions with a best-fit corresponding to 
maximal mixing. This remains true also when the K2K data is included or separately analyzed 
(e.g., ^ni)- Though exactly maximal mixing and zero Ue3 would hint to some underlying 
symmetry in the lepton sector, one can not expect radiative corrections to allow these extreme 
values to persist down to low energy ^28^. Thus, one expects non-extreme values for 6*13 and 
623- See, e.g., [231221 for ways to generate deviations from the bimaximal mixing scheme. 
All in all, the three observations (i) to (iii) together with the mixing matrices (|8I9|) lead us to 
parametrize three elements of the mixing matrix as 

U,3 = AX^, (12) 
f^M3 = y^(l-5A-)e*^ • 

For A = one would have the bimaximal scheme from Eq. The two Majorana phases 

are left out for the moment. Unitarity of t/pMNS suffices to calculate the remaining elements. 
The parameters A and B are numbers of order one. The 623 ^ 7r/2 — 6*23 ambiguity reflects in 
a sign ambiguity of B. The power of A in the expressions for 11^3 {U^s) can be adjusted when 
more stringent limits (more precision data) are available. 

We can take the best-fit value from Eq. (jHI) to calculate A ~ 0.22, which is remarkably similar 
to the Wolfenstein parameter or the sine of the Cabibbo angle. The maximal allowed value of 
l^esP = 0.07 corresponds to Ue3 = AX with A ~ 1.2. At 1 (3)cr, the range of A lies in 

A ~ (0.08)0.18 -0.28(0.35) . (13) 

For the best-fit points of the many available analyzes [H], which lie in the range between 0.41 
and 0.44 for tan^ 6^12, A is between 0.24 and 0.22. 

If indeed A ^ 0.22, then for m = 1 it must hold that B ^ 0.91 in order to fulfill the requirement 
sin^ 2^23 ~ 0.85. In the following we shall work with the "best-fit" value of A = 0.22. If the 
limit on |t/e3p goes below ~ 10~^ one should take the power n = 2 in Eq. (fT^ in order to keep 
A of order one. If |f/e3p ~ 10~^, then n = 3 is advantageous to choose. 

Analogously, since typically (see below) sin^ 26'23 — 1 — 45^ A^"* one should for values larger 
than sin^ 26*23 — 0.95 (0.99) use the power m = 2 (m = 3) in Eq. (jl2j) . Values of m = 4 

^Deviation from maximal solar neutrino mixing has also been analyzed in terms of the parameter e = 
1 — 2 sin^ 612 [317, which roughly corresponds to A: e ~ A + O(A^). 



6 



would be required if a precision in sin^ 2^23 of order 10"'* was present, which seems improbable 
unless a neutrino factory will be operative. In terms of tan^^23, which in the future will be 
more appropriate to use, one will find that tan^ 6^23 ~ 1 — 4 5 A™. Thus, for tan^6'23 ^ 0.7 (or 
^ 1.3) one should take m = 2, while for tan^ 6*23 ^ 0.9 (or ^ 1.1) the value m = 3 is more useful. 

We shall now consider several different cases for the powers of A in Eq. (jl2j) . The considerations 
from this Section indicate that current data and the precision of future experiments on ^23 and 
^13 limit the realistic values of m and n between 1 and 3. 



3 The Mixing Matrix 

3.1 Case m = n = 1 

In this case we have Ues = AX and U^s = -y/|(l ^ B X) e^^ . It corresponds to rather large 
deviations from the extreme bimaximal values. One "predicts" f/es very close to its current 
limit and also sin^ 2^23 is on the edge of its 3cr range. We identify 

tan2 ^12 ^ 1 - 4 A + 2 (5 + ^2) A2 + ^(A^) , 
sin^ 2^23 ^ 1 -452 A2 + 0(A3) , tan^ ^23 ^ 1 - 4 5 A + 2 (A^ + 5 5^) A2+ ^(A^) (14) 

sin^ = ^2 A2 . 

The fact that sin^ 2^^23 is blind to the ^23 ^ 7r/2 — ^^23 ambiguity is reflected in the fact that B 
appears quadratically in the last expression. The form of ?7pmns is rather lengthy, and we shall 
give it therefore only to order A: 



f^PMNS 



I yj(l + A) yj(l-A) ^ 

- {I- B - Ae'^)\) \{l + {l + B- Ae'^)X) J\{l-B\)e'^ 



+ 0{X^ 



\ \{l-{l + B + Ae'^)\) -\{l + {l-B + Ae'^)X) ^\{l + B\)e'^ ) 

(15) 

The precise form to a given order of A is easily obtained by using the unitarity of the mixing 
matrix. The corrections quadratic in A are functions of A and B except for Uei-, which receives 
only corrections depending on A. Important to note is that the corrections from A are respon- 
sible for the deviations from the values ±1/2 of the entries in the lower left 12 block. Finally, 
the invariant measure oi CP violation in neutrino oscillations is 

Jcp = ^ (1 - (2 + + 2 B^) X^) sin 5 + O(A^) . (16) 

Noting that AX ~ ^13, the corrections stemming from ^13 and non-maximal ^12,23 — as indi- 
cated in Eq. (fTIHl — are easily identified. The larger the deviations from maximal 6'i2,23, i-e., 
the larger A and B, the smaller becomes Jcp- The "prediction" is that CP violating effects 
are up to ~ 5%. Note however that actual experiments searching for leptonic CP violation will 
not just measure Jcp, see Section EIH 
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3.2 Case m = 1 and n = 2 

For these values we have 11^3 = AX^ and U, 



fi3 



Ul- BX)e'^. We can identify 



tan2 6l 



12 



1 - 4A + 10A2 + C»(A3 



sin' 2^23 -l-4:B'X' + 0{y) , tan' ^23 ^ 1 - 4 5 A + 10 5' A' + O(A^) 



(17) 



sin^ = A^X^. 



The "predictions" are \U, 



e3| 



10 ^ and atmospheric mixing very close to the end of its allowed 



3a range. The mixing matrix f/pMNS reads 



f^PMNS 



i(l + A) 



Hi-A) 



AA^ 



-|(l-(l-fi)A) i(l + (l + 5)A) JUl-B\)e 



y 1(1 -(1 + 5) A) -i(l + (l-fi)A) ^lil + BX)e'^ J 



It is obtained by removing the term A e^^ from the PMNS matrix in the case of m = 77. = 1 as 
given in Eq. (fT3j) . The corrections of order A^ for the lower left 2 by 2 submatrix are functions of 
A and B. They are constant for f/gi and only depending on B for Effects of CP violation 
are proportional to A^, 



Jcp 

and not more than a few %. 



AX' 



{1 - 2 {1 + B^) X^) sm6 + 0{X^ 



(19) 



3.3 Case m = 2 and n = 1 

Now our parameters read Ue3 = AX and = ^/|(1 — B A^) e*"^. The mixing angles are 

tan^ ^12 ~ 1 - 4 A + 2 (5 + ^2) A^ + ^(A^) , 
sin2 2023-l-(^^-2 5)2A^ + O(A5) , tan^ ^23 - 1 + 2 (A^ - 2 5) A^ + ^(A^) , (20) 

sin^ ^13 = A^X' . 

Thus, Ue3 is close to its current limit and the deviation from sin^ 2^23 = 1 is not more than a 
few %. The mixing matrix f/pMNS is given by 



t^PMNS 



Ul + X) 



Hi-A) 



AX 

-i (1- (l-v4e^'^)A) i(l + (l-Ae^^)A) (1 - B X^)e'^ 

|(l-(l + Ae^^)A) -i(l + (l + Ae^^)A) ^^{l + B X')e'^ j 
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+ 0{X^) (21) 



which is obtained from Eq. (jl5j) by removing B from the lower left 2 by 2 submatrix. The 
quadratic corrections are functions of A and B except for Uei, which only depends on A. The 
rephasing invariant CP violation measure is 



J. 



CP 



AX 



(1- {2 + A^)X^) sm6 + 0{X' 



being rather sizable but not exceeding 5%. 



(22) 



3.4 Case m = n = 2 

Now it holds U^^ = AX^ and U, 



B A^) e*^, yielding the mixing angles 



tan^ ^ 1 - 4 A + 10 + ^(A^) , 
sin^ 2^23 ~ 1 - 4 52 A^ + 0{X^) , tan^ ^23 ^ 1 - 4 S A^ + ^(A^) 

sin^ = ^2 A^ . 



(23) 



The deviation from sin^ 2^^23 = 1 is not more than a few percent and |?7e3| is on the level of 
10~^. The mixing matrix f/pMNS is given by 



f^PMNS 



i(l + A-A2) 



^(1-A) 



AX'' 



i(l- A+ (S + Ae*^) A2) 1 (i + A- (l-S + Ae*^)A2^ 



{l-BX')e 



\ i(l-A-(5 + Ae^'5)A2) -1 (i + A-(l + 5-Ae^'5)A2) ^ ^(i + B X'^) e'^ J 



+0{X^ 



(24) 



It is seen that for the lower left 2 by 2 submatrix the linear corrections to the "bimaximal" 
values ±1/2 are constant and the quadratic ones are functions of A and B. The rephasing 
invariant CP violation measure is given by 



J, 



AX' 



CP 



(1-2A2) sm6 + 0{X^) 



(25) 



again on the level of a few percent. 



3.5 The "Wolfenstein Case" m = 2 and n = 3 

How could one not be tempted to put the third power of the expansion parameter in the Ues 
and the second power in the ?7^3 element. This would resemble the Wolfenstein parametrization 

Eq. dH). In this case, i.e., Ues = AX^ and U^3 = •y/|(l — B A^) e*"^, we have 

tan2 ^ 1 - 4 A + 10 A2 + ^(A^) , 
sin^ 2^23 ^ 1-452 A^ + 0(A5) , tan2e23 ^ 1-45A2 + C(A3) , (26) 

sin^ = A''X^. 
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The "prediction" for sin^ 2^23 — 1 is not more than a few percent and iC/esI^ is on the level of 
10~^. For the PMNS matrix one has 

C^PMNS — 

'i(l + A-A2 + A3) yj(l_A) AX^ 

{l-X + BX'^-{B-Ae'^)X^) ^ {l + X - {1 - B) X'^ + {1 + B - Ae'^) X^) ^{l-BX^)e'^ 

\ \ {l-X-BX^ + {B- Ae^^)X^) -\ {l + X- [1 + B)X^ + {I- B + Ae'^)X^) ^{l + BX^)e'^ ) 

+0{X^) . 

(27) 

In contrast to the Wolfcnstcin parametrization, however, A appears linearly in e.g. the Ufj,2 
element. Also, CP violation is proportional to the third power of A, 

A A3 

Jcp ~ ^ (1 - 2 A') sin 5 + C»(A^) , (28) 

not exceeding one percent. 

3.6 The remaining cases 

The remaining interesting cases are m = 3 with n — 2 and also m — n — The former case 
— that is = (1 - 5 A^) e'^ and Uez = A A^ — yields 

tan2^i2 ~ 1-4A+10A2 + C>(A3) , 
sin2 2e23- l-4S2A6 + e)(A^) , tan^ ^23 ^ 1 - 4 S A^ + O(A^) , (29) 

sin^ ^A^\\ 

with \Uez\^ ~ 10~3 and sin^2^23 7^ 1 with a precision of less than 1%. For the PMNS matrix 
holds: 

t^PMNS — 

( Y^(l + A- A2 + A3) \ 

-\ {l-X + Ae'^X^ + {B + Ae^^)X^) \ (l + A - (1 + Ae*'^) A^ + (1 + 5 + Ae^^) A^) ^{l-BX^)e^^ 

^ \ [l-X- Ae'^ X^ -{B + Ae^)X^) -\ {l + X- {I - Ae'^) X^ + {I - B - Ae'^) X^) ^{l + BX^)^^ j 

+C»(A^) . 

(30) 

As for the case m = n = 2 one obtains: 

A A2 

Jcp ~ ^ (1 - 2 A^) sin S + 0{X') , (31) 
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again at most a few percent. 

If m = n = 3, i.e., ?7^3 = — B X^) e''^ and Ues = AX^, then it holds for the mixing 

parameters: 

tan^ ~ 1 - 4 A + 10 + 0{X^) , 
sin^ 2^23 ^ 1-452 A6 + 0(A^) , tan2 023 ^ 1-45A3 + C(A^) , (32) 

h3 



sin^^i, = ^2 A^ 



i.e., except for sin^ ^13 to the given order identical for the m = 3 and n = 2 case above. The 
parameters A and B appear in the mixing matrix only a third order in A: 

f^PMNS — 

/ yi(l + A- A2 + A3) ^(i-x) AX^ \ 

{1- \+ {B + Ae'^)X^) i (l + A-A2 + (l + 5-^e^^)A3) ^(i-BX^)e'^ 

\ ^ {l-X-{B + Ae'^)X^) {l + X-X'^ + il-B + Ae'^)X^) + B X^) e^^ J 

+0{X^) . 

(33) 

There are no quadratic terms in the elements U^i and Uri- Finally, CP violation is governed 
by 

Jcp ~ ^ (1 - 2 A^) sin 5 + 0(A6) (34) 

being below one percent. 

The remaining cases m = 3 and n = 1 (m = 1 and n = 3) are obtained from the cases m = 2 
and n = 1 (m = 1 and n = 2) by setting S = (A = 0) in the relevant expressions for the 
mixing parameters and Jcp- 



3.7 Speculations 

One can speculate about the origin of the corrections induced by the A terms. It is possible to 
imagine, e.g., that the bimaximal mixing scheme from Eq. (|lip stems from the diagonalization 
of the neutrino mass matrix (this is possible, e.g., when a — — symmetry is present 
PTj ) and any corrections are implied by the unitary matrix Ui that diagonalizes the charged 
lepton mass matrix [27\. Recall that in a basis in which the charged lepton mass matrix is not 
diagonal the PMNS matrix is given by Uj U, where U diagonalizes the neutrino mass matrix in 
that basis. If we define the matrix Ux, which induces the correction to the bimaximal scheme, 
we may write ?7pmns = f^Af^pMNS; where t^pMNs given in Eq. (fTTjl ^. Then one can simply 

^This is similar to the strategy in where however a different mixing matrix to start with was used. 
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solve for Ux. Taking for definiteness the example m = 3 and n = 2, one finds 



/ 



1 - AV2 





v 



2 + X ~^ T J 



+ 0{X^) . (35) 



It is seen that 23 and 32 entries are in general of order one but reduce to order A^ for CP 
conservation. Those entries can also be of order A, however only for for the cases m = n = 1 
and m = 1 with n = 2. Thus, if CP is conserved and atmospheric neutrino mixing is very 
close to maximal, the matrix Ux takes the unit matrix as the dominant form with corrections 
of order A. The typical "CKM-structure" with very small A^ terms is however not necessary. 

4 The Majorana Mass matrix 
4.1 Basics 

Up to know our analysis assumed only the neutrino oscillation explanation of the experimental 
data. Now we assume in addition that neutrinos are Majorana particles, which is, e.g., a 
prediction of the see~saw mechanism [S21- Thus, the neutrino mass matrix m,^ in the basis in 
which the charged lepton mass matrix is diagonal is given by: 



Here m'^^''^^ is a diagonal matrix containing the masses mi_2,3 of the three massive Majorana 
neutrinos. Immediate consequence of the Majorana nature of the neutrinos is the presence of 
two Majorana phases a and (3 to which neutrino oscillations are insensitive Information 
about these phases can be obtained by studying processes in which the total lepton charge L 
changes by two units, e.g., neutrinoless double beta decay, — >• tt" + /i"*" + /i"*", etc. Realisti- 
cally, only neutrinoless double beta decay can expected to be measured |34l- The decay width 
of this process is sensitive to the ee element of m,y. 

An interesting observation is that the ratio of typical best-fit values of the mass squared dif- 
ferences corresponds roughly to the expansion parameter A: 



We took for Am\ the best-fit point of the preliminary new analysis of the SuperKamiokande 
collaboration [SO]- Using for Am\ the 90 % CL. analysis from [20^, which is (1.3 . . . 3.1) ■ 10~^ 
eV^, with the 90 % CL. range of AttIq from j^I, we find that R lies between 0.13 and 0.28. 
This corresponds to a good precision to the 3a range of A as given in Eq. ()13|) . In the following 
we shall assume that i? ~ A and study the resulting structure of the neutrino mass matrix. 




PMNS 



(36) 




(37) 
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The results do not change much unless Am\, Amg and tan^ 612 are on the very edges of their 
allowed ranges. Before we perform this analysis, it is useful to study the mass matrix again in 
the limit of exact bimaximal mixing. In the following, we will neglect the CP violating phases, 
see, e.g., for an analysis of the structure of uiy in case of complex entries. Using Eqs. Q 
and the mass matrix reads 



nit, 



( A 



\ 



B 

2 



-B 


\ 










2 








Z) + ^ 

2 


/ 



(38) 



where 



A 



nil + ni2 



B 



ni2 — nil 
2V2 



D 



ni3 
2 



(39) 



As mentioned, there are three extreme cases for the mass hierarchies, the normal hierarchy (NH) 



with a/At 



m^., 3> mi 



0, the inverse hierarchy (IH) with a/ Am 



im^ = ms 3> 7712 

m2 — mi ^ 7713 ~ and quasi-degenerate neutrinos (QD) with m^ = m^ ~ m2 — mi. 
Depending on the relative signs of the mass states, several extreme forms of the mass matrix 
result. In case of NH, one finds for mi 2 = 0: 



/ \ 

■ 1 1 

[■ ■ l) 



(40) 



i.e., a mass matrix with a leading /ir block. Regarding IH, the third mass m^ can safely be 
neglected. The form of mi, then depends on the relative sign of the two mass states mi and 
7712: 

' / 1 \ 



m 



1 _i 

2 2 



v 



/o 1 

I V • • 



V2 V2 \ 





same sign 



opposite sign 



(41) 



/ 
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For the QD spectrum one finds 



rriiy = mo < 



v 

/o 

v ■ 

/o 
V ■ 



_ 


\ 


1 







1 / 


1 


1 


V2 


V2 


1 


1 


o 
Z 






1 


1 


2 
1 






1 


1 


2 


2 




1 




2 





^ 





-1 



sign(mij = sign(m2j = sign^msj 



sign(mij = — sign(m2j = sign^msj 



(42) 



sign(mij = — sign(m2j = — sign(m3j 



I v 



sign(mij = sign(m2j = — sign^msj 



/ 



We can expect that in our parametrization the parameter A will appear in the neutrino mass 
matrix at least linearly in order to correct the extreme values 0, ±1 and ±1/2. 

4.2 Normal hierarchy 

In case of the normal hierarchy we have 



Am\ , 1712 



Am^ A and mi = D 



Ami 



where / > 



(43) 



The expression for mi with D = 0{1) expresses our lacking knowledge of it. A similar Ansatz 
for the structure of m^, in case of a normal hierarchical mass scheme has been made in [TT] . For 
m = n = I = 1 and all mass states positive the mass matrix looks like 



( A 



V 



1+2 A 



A 



2A-1 
V2 



A 



1 + 



2B\ A 



(44) 



l + (i+2S) A / 



neglecting terms of order (9(A^). The characteristic "leading /ir block" structure of m^ from 
Eq. ()4()j) is identified. Corrections at order A depend on A in the e-row of m^ and on B for the 
IHX and rr entries. Higher powers of A in mi will — to order A^ — lead to the disappearance 
of D in the formula for m^. Setting ?7e3 = A}? leads to a mass matrix in which to order A^ the 
parameter A does not appear in the [it submatrix as well as in the ee entry. It is obtained by 

removing A from the last equation in the indicated entries. The matrix for = w ^ (1 — i? A^) 
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but Ue3 = AX is given by Eq. (j44j) by removing B from the first row and from the hnear 
terms of the fifi and rr entries. The "Wolfenstein hke" parametrization with Ues = AX^ and 
1 



/x3 



1 — B X ) together with mi 
/ A - 2 



X 

V2 



D a/ Am\ X^ leads to the particularly simple form 



1 + | + (1-2 5)A2 



A 

l-A^ 



v 



+ 0{X^ 



(45) 



1 + f + (l + 2fi)A2 / 



The 023 ^ tt/2 — 023 ambiguity, which translates into a sign ambiguity of B, is seen to have 
origin in the size of the /i/i and rr entries. E.g., for all masses positive and 623 > 7r/4 the 
rr entry is larger. The structure of the mass matrix does not depend on the exactness of 
the relation i? = A or the relative signs of the mass states. When sign(mi) = — sign(m2) = 
— sign(m3), then the mass matrix looks as above. For sign(mi) = — sign(m2) = sign(m3) and 
sign(mi) = sign(m2) = — sign(m3) one has to replace^ A with —A, B with ~B and the 1 in the 



Hi 



B X^) than we obtain a mass matrix in 



/IT block with —1. If we further choose U^3 = 

which up to order A^ the parameters A, B and D do not appear at all. It is obtained by setting 
in the last equation B = 0. 



4.3 Inverse hierarchy 



In this case we have 



m2 



Am 



A ' 



mi 



Am\ 



1 - AV2) and mg = D 



Ami ^2+' 



where / > . 



(46) 



The dependence on the power of A in 1713 is almost vanishing. The form of iriy depends strongly 
on the signs of the masses mi and m2. For identical relative signs between mi and m2, the 
e/i and er entries are suppressed by A or A^, depending on the powers of A in or Ue3. If 



e.g., m = 2 and n = 1 or m 



n 



1, then the entries are of order A. For all other cases 



under consideration, these terms are of order A^. The remaining independent entries of m,^ are 
order one. If mi and m2 have opposite relative signs, then me^ and mg^ are order one and the 
remaining entries of m,^ are linear in A, independent of m and n. One finds for m = n = 1 that 
for same signs of mi and m2 



/ 1 



m. 



Am\ 



l + BX 



4 + i?A 
l-BX J 



+ 0{X^) 



(47) 



while for opposite signs 



m,. 



Ami 



/ 2A 



V 



1 

V2 



1 

V2 



(A-l)A 



A 



0{X' 



(4J 



-{A + 1)X J 



^The convention here and in the following will be such that the sign of the ee entry is positive. 
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The parameter A appears at least linearly to correct the extreme "bimaximal" mass matrices 
from Eq. ()41|) . For same signs, A appears in the e-row and B in the /ir sector, whereas for 
opposite signs it is vice versa. Taking as another example again the "Wolfenstein-like" case 
m = 2 and n = 3 one finds to order A^: 



( l-U' 



1 \2 
i72 ^ 



( 2A 



V 



5 \2 



1 I 9-4 B 1,2 



-{B + \)\' ) 

J 9+4 B \2 \ 

V2 4v^ ^ ^ 



A-fA^ 
-A + lA^ 



same signs 



opposite signs 



(49) 



As usual, B will not show up for higher orders of A in [7^3, leading to a mass matrix that is at 
order A^ only a function of A. It is again obtained by setting in the last equation i? = 0. 



4.4 Quasi— degenerate Neutrinos 



For quasi-degenerate neutrinos, i.e., m| ~ m\ ~ m\ ^ Am^, there is another small quantity 
introduced, namely the ratio of the common mass scale mo with /S.m\. For simplicity we work 
with the normal mass ordering. In this case we can express the three mass states as 



= iriQ , m2 = a mo and mi = b mg , 
where a = 1 — v , b = 1 — r] {1 + A^) 



and rj 



(50) 



2 mi 



The common mass scale is denoted by mo. These expressions for the masses are valid to order 
Tj. Since the spectrum is quasi-degenerate for mo ^ 0.2 eV, we can estimate rj ^ 0.04 eV, 
therefore \^ > rj. 

First, we take the case that all mass states have the same relative sign. For m = n = 1 we find: 



m. 



mo 



/ 1-r/ 



v 



l-l-BriX 



Taking the case m = 2 and n = 3 we have 



fir/A / 



(51) 



/ 1 



m,. 



mo 



ri2?_ 

1+4 B 



2^2 

2 _|_ '!l2^ 

R _|_ 4g-l 
2 4 



+ 0{r]X^ 



(52) 
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It is seen that for rj = the mass matrix is proportional to the unit matrix, irrespective of m 
and n. The corrections to the zero values of the extreme bimaximal form from Eq. (jl^ are 
very small. 

If sign(mi) = — sign(m2) = sign(m3), then the dependence on i] is not so important. Neglecting 
7] with respect to terms of order 1, we find for m = n = 1: 



/2A -:^ + 



1 I A-B 



im,. 



mo 



A 



1 _|_ A-B 



\ 



l + iA-B-l)\ i + A 

i - (A + 1 - 5) A y 



(53) 



and for m = 2, n = 3: 

/ 2A - A2 

= mo 



[m 



l-T) , 2~B \2 



1-ri _ B+2 \2 
V2 V2 ^ 

I + A - I A^ 



V 



■ i-A+(|-5)A^ 

i-A+(i + i?)AV 
For the case sign(mi) = — sign(m2) = — sign(m3) one finds for m = n = 1 



(54) 



mo 



/ 2A 



v 



1 A+B 



X 



V2 V2 

1 + {A + B-1)X 



J A+B 

V2 V2 

-k + x 



+ 0{r],X^' 



(55) 



{A + B + 1)X J 



and for m = 2, n = 3: 

/2A-A2 -i^ + ^A^ 



(m^)+__ = mo 



V2 V2 

|-A+(| + 5) A^ 



V 



B±2 \2 \ 

-\-X+[\-B) X- I 



+ 0(r/A,A3) 



(56) 



These two last cases look very similar. Finally, the situation for sign(mi) = sign(m2) 
— sign(m3) looks simpler: e.g., when m = n = 1: 



/1-2A2A2 -^AX + V2ABX'^ -^plAX-^flABX^ \ 



m. 



mo 



2BX-B^ A^ 



V 



-1 + {A^ + 2B^)X^ 
-25A+ (2^2 + ^2) A2 / 



o{v,x') . 



(57) 



The corrections to the entry 1 (0) that is present in the extreme bimaximal form from Eq. 
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are at least quadratic (linear). For m = 2 and n = 3 it holds: 





(m^)++_ = mo 



2+(i?(2-r/)-|) 



(5(2-r^) + 2) ; 



+ 0(A3) 



V 



2 

2 



(58) 



4.5 Summary for the Mass Matrix 

Looking at the cases considered in the last Subsections, the following summarizing statements 
can be made: 

• Roughly, for |?7e3| ~ 0.01 and sin^ 2^23 ~ 0.9, corrections to the extreme forms of the 
mass matrices in Eqs. ()40I41I42|) are linear in A. When |?7e3| ^ 10~^ and sin^ 26*23 ~ 0.95, 
the corrections become quadratic. 

• For the normal hierarchy, corrections to the exact bimaximal form are at least order A. 
To lowest order, the parameter A appears in the e row of m^, and B in the /i/i and rr 
elements. The number D, parameterizing the unknown lightest mass state, appears in all 
entries. There is basically no dependence on the relative signs of the mass states. 

• For the inverse hierarchy, the dependence on D is highly suppressed. For identical signs of 
the two heaviest mass states, the correction to the ee entry, whose extreme value in case 
of bimaximal mixing is (in units of Am\) 1, is at least order A^. The remaining elements 
receive at least linear corrections. A appears at leading order in the e-row and B in the 
fiT sector. Opposite signs of the two leading mass states lead to linear corrections to the 
entries and the appearance of A in the fir sector and B in the e-row. 

• In case of a quasi-degenerate spectrum and identical signs of the masses, the corrections 
to the unit matrix are at least quadratical. The cases sign(mi) = — sign(m2) = sign(m3) 
and sign(mi) = — sign(m2) = — sign(m3) look very similar. For sign(mi) = sign(m2) = 
— sign(m3) there are only quadratic corrections to the ee and fir entries. 

• The /ir entry is special in our parametrization since the parameters A and B do typically 
only appear there for rather large deviations from zero Ues and from maximal atmospheric 
neutrino mixing. 

• If we would consider the inverse ordering in the QD mass spectrum, one has to change 
Eq. (jHn|) to 7712 = fTT'O, fTT'i = cL7nQ and 771^ = 6 mo. Only simple sign changes for A and/or 
B in some elements of m^, would be the result. The main difference would be for the case 
sign (mi) = sign(m2) = sign(m3), where the corrections to the e/x and er elements are 
now order 7] and not just 7] X or 7] X"^. 
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5 Applications 



It is surely useful to study formulae which are obtained by expansions in small parameters or by 
certain simplifications within our parametrization. We shall perform this analysis now for the 
effective mass as measurable in neutrinoless double beta decay and the oscillation probabilities 
for long-baseline neutrino oscillations. 



5.1 Neutrinoless Double Beta Decay 



We shall analyze now within our parametrization the form of the ee element of m^, which 
is denoted by (m). In a given mass scheme or hierarchy one can considerably simplify the 
expression for (m) We first note that since 



the results are independent on the power of A in 



(59) 



5.1.1 Normal Hierarchy 



With the help of Eqs. ()43I59|) we can evaluate the effective mass in case of the normal hierarchy. 
We find for Ues = AX and 1713 = A^: 



(m) 



A-2A^ 1 



D 



C2a — A? C2(a-I3) 



+ C(A^ 



(60) 



where C2a = cos 2a and so on. Terms of order A^ depend on A, D and the two Majorana phases. 
Choosing = D or higher powers of A leads to the disappearance of D in the formula. 
For m3 = D X^ and higher orders of A in Ug^, i.e., n > 2, it holds 



(m) 



Am 



X-2X^ (l-^^C2a) +^ (4(1 + 2Dc2.)+D2(i_c4„)) 



+ 0(A^) . (61) 



The formulae for 777,3 = D X^ correspond to setting = in this last equation. Roughly, we 
can estimate the effective mass in the normal hierarchy as 



{m) 



X 



Ami ^ 0.005 eV . 



(62) 



5.1.2 Inverse Hierarchy 



From Eqs. ()46p59j) one sees that in the expression for (m) the dependence on j3 practically 
vanishes. The result for Ues = AX is 



{m) 



Ami W 4 + 7 (7 - 4 - (9 + 4 A^)c2a) X^ + O(A') 



(63) 
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Higher powers of A in f/gs lead to the disappearance of A in this equation. The maximal and 
minimal values are obtained when a takes the values and 71/2, respectively. Thus, 



2 yX^A < (m) < - (^A' + A^^ , 

up to corrections of 0{X^). For no extreme values of a, the scale of the effective mass is 



(64) 



(m) ~ y Ami ^ 0.05 eV . (65) 

Comparing this with the value of (m) in the normal hierarchy in Eq. (j62j) one sees that the 
expansion parameter A shows up as the ratio of the typical values of (m) in the inverted and 
normal hierarchy. 

It is known that extraction of information from a measurement of 0^(3(3 suffers from a large 
uncertainty stemming from the calculation of the nuclear matrix elements. This uncertainty 
is a number of order one jHZI. It is therefore an important question to ask and an even more 
important one to answer whether future Oi>f](3 experiments can distinguish pHj l^ between the 
normal and inverted mass hierarchy. Let us parametrize the nuclear matrix element uncertainty 
with a factor ^ as done in [211 • In order to distinguish the normal from the inverted hierarchy 
it must hold that the maximal value of (m) in the normal hierarchy times the uncertainty ^ 
has to be smaller than the minimal value of (m) in the inverted hierarchy. Therefore, choosing 
Ues = AX and small 1713 we find from Eqs. ()60I64|) 

e<4 (1 + 2A(1±A2)) +0(A=^) . (66) 

Needless to say, A vanishes for smaller values of Ue3- If that is the case, then ^ 6, which is 
a very realistic number. Thus, with our expansion parameter A ~ 0.2 and |f/e3p ~ 0.01 it is 
easily possible to distinguish between the normal and inverted mass hierarchy. 

5.1.3 Quasi— degenerate Neutrinos 

The formulae for the mass states are given in Eq. (j43j) . Ignoring rj and taking 11^3 = AX one 
finds 

^ ^ \/cl + (2 - - (2 + A^) C2a + 2 C„ C„_2/3) A^ + 0(A=^) . (67) 

mo V 

Interesting are the cases which correspond to CP conservation, and which are obtained by 
setting a, /? to 7r/2, vr. They read up to OirjX^, X^): 

1 — rj a = /3 = ^ sign(mi) = sign(m2) = sign(m3) 

2A 2a = (3 = ^ sign(mi) = — sign(m2) = sign(m3) 

2 A a = (3 = 7r/2 ^ sign(mi) = — sign(m2) = — sign(m3) 

^ 1 — ?7 — 2A^X^ a = 2f3 = TT sign(mi) = sign(m2) = — sign(m3) 



(m) 
mo 



(6J 
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As usual, for = AX'^ and above, the dependence on A (and thus also on (3) drops and 
appears only to order A'*. Noting that the minimal value of (m) is 2Amo, we can investigate if 
future experiments can distinguish between the normal and quasi-degenerate mass hierarchy 
jSni- In analogy to the discussion leading to Eq. (jUBj) . it follows from Eqs. ()6UI68p 

e< (1 + 2A(1±A2)) , (69) 

which will be easily possible. It is a bit more tricky to distinguish between the quasi-degenerate 
and the inverted hierarchies. The requirement for ^ is 

e<yjA(l+(A^ + l)A^) . (70) 

which is suppressed roughly by a factor A with respect to the limit on ^ in order to distinguish 
the normal and quasi-degenerate mass hierarchy. Also in this aspect the parameter A shows 
up as a scaling factor. 



5.2 Long Baseline Oscillation Experiments 

There is another field of neutrino physics in which expansion in small parameters gives insight 
in the physics involved and which is therefore useful to study within our parameterization. 
These are the oscillation probabilities for long-baseline experiments . The determination of 
some currently unknown neutrino parameters, namely f/es, the sign of Am\ and the Dirac-like 
CP violating phase are the purpose of such experiments. There are helpful expansions of the 
relevant oscillation probabilities in vacuum [IHl- Here, we do not consider matter effects since 
they will not change our conclusions. Let us first comment on CP violation. Using Eq. (^, 
one finds for the difference of the oscillation probabilities: 

AP = P{Ve ^ V^.) - P{Ve ^ V^) = 

(71) 

which, using Am\2 = Amg^ — Am\]^, can easily be shown to vanish for two masses being equal. 
The invariant Jqp was defined in Eq. (fTUI) . Since Am\^/ Am\^ = ±A^, where the "+" in case 
of normal ordering and the "— " for inverse ordering, we can expand the last equation: 

AP ~ ±4 JcP X' sin2 + OiX^) . (72) 

Thus, the CP violating effects in realistic experiments are suppressed by another two orders of 
A in addition to the suppression present in Jcp- If n is the power of A in Ues, then the total 
suppression is A^"*"". 

One can also consider the bare oscillation probability for the "golden channel" , which is given 
by i^e ^fj. oscillations. Using the form of P{ve ~^ ^pi) as given, e.g., in [22], one finds for the 
oscillation probability in case of m = = 1: 

P{ve ^ z/^) ~ 2 A^ sin^ A32 A^ - 2 A sin^ A32 {2AB + cos((5 T A32)) A^ + C»(A^) . (73) 
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Here, the "— " sign is for neutrinos and the + for antineutrinos. We defined A32 = (rrig — 
rn^L/AE. The first term proportional to is the term that probes Ue^ whereas the second 
term proportional to A^ is the one probing the CP phase 5. 



As another example, assume n = m = 2. Then, the terms probing Ug^ and CP violation will 
both be proportional to A"^: 

P(i/, ^ u,) ^ ^^^^ {1 + 4A'-4A cos(5 T A32)) + 0{\') . (74) 

The parameter B only appears at order A^, since there are no terms of order A^. For m = 2 
and n = 3 it holds: 

P{u, ^ u,) ~ ^Jll^ X^-2A sin^ A32 cos(5 T A32) A^ + O(A^) • (75) 

A characteristic combination of the oscillation parameters that appears in the relevant proba- 
bilities is Am2i/Am|^ sin 2^12 [10] • Neglecting terms of order A^, we find for this parameter in 
our parametrization that 

;i -2A2 + 2(1 + A2)A3) forn=l 

Arn2 

sin 2^12 c^±\^ { (1 - 2 A2 + 2 A3) forn = 2 , (76) 



Amii 



1 - 2 A^ + 2 A^) forn = 3 



where again the "+" in case of normal ordering and the "— " for inverse ordering. The difference 
between the cases n = 2 and n = 3 appears only at seventh order in A. The characteristic 
parameter is therefore to order A"^ independent on the precise form of the parametrization. 



6 Conclusions 

The zeroth order approximation for neutrino mixing can be the bimaximal scheme with two 
maximal and one zero angle in the mixing matrix. It can be used as a reference matrix, whose 
corrections can be described in a similar manner as the Wolfenstein parametrization describes 
corrections to the unit matrix. Indeed, at least one of the angles in neutrino mixing is different 
from the extreme value corresponding to bimaximal mixing, namely the angle describing solar 
neutrino oscillations. To take this into account, a flexible parametrization of the neutrino mixing 
matrix was proposed in which the expansion parameter A ~ 0.2 is introduced to quantify this 
deviation from maximal mixing of solar neutrinos. It can also be used to quantify the possible 
deviation from zero f/gs and maximal mixing of atmospheric neutrinos. The power of A to 
usefully describe these two latter aspects can be adjusted to future data. Depending on the 
power of A, rather simple forms of the PMNS matrix are obtained, where the deviations from 
the "bimaximal" values 0, ±1/2 and ±l/\/2 are implied by A. If ?7^3 and t/gs are close to 
their maximally allowed values, A appears at first order in all elements of f/pMNS- For values of 
If/gsl ~ 10~3 and sin^26'23 ^ 0.95, the corrections become quadratic. The invariant measure for 
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leptonic CP violation is proportional to A", where n is the power of A in Ues- One can interpret 
these corrections to the exact bimaximal mixing scheme in the same way as corrections to the 
unit matrix lead to the CKM matrix for the quark sector. Observing further that the ratio 
of the mass squared differences as measured in experiments is roughly A^, allows to study the 
form of the Majorana neutrino mass matrix m^. Also here, the corrections to the extreme 
forms of nil, in case of bimaximal mixing and extreme hierarchies are linear or quadratic in A, 
depending on the precise values of Ues, Uf^s or the value of the smallest mass state. The ee 
element of m„ can be measured in experiments probing neutrinoless double beta decay. Here, 
A appears as the scale factor of the typical values of (m) in the normal and inverted hierarchy. 
It also influences the maximal value of the uncertainty in the calculations of the nuclear matrix 
elements allowed to distinguish the normal, inverted or quasi-degenerated mass hierarchies. We 
furthermore commented on how our parametrization applies to realistic long-baseline oscillation 
experiments. Simple forms of the relevant oscillation probabilities are obtained. In particular, 
due to the small ratio of the two independent mass squared differences, effects of CP violation 
are suppressed by another two orders of A. 
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